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We derive tlie Ornstein-Zernilie equation for molecular crystals of axially symmetric particles and 
apply the Percus-Yevick approximation to this system. The one-particle orientational distribution 
function p^^\VL) has a nontrivial dependence on the orientation Q,, in contrast to a liquid, and is 
needed as an input. Despite some differences, the Ornstein-Zernike equation for molecular crystals 
has a similar structure as for liquids. We solve both equations numerically for hard ellipsoids of 
revolution on a simple cubic lattice. Compared to molecular liquids, the orientational correlators 
in direct and reciprocal space exhibit less structure. However, depending on the lengths a and b 
of the rotation axis and the perpendicular axes of the ellipsoids, respectively, different behavior is 
found. For oblate and prolate ellipsoids with 6 > 0.35 (in units of the lattice constant), damped 
oscillations in distinct directions of direct space occur for some of the orientational correlators. They 
manifest themselves in some of the correlators in reciprocal space as a maximum at the Brillouin 
zone edge, accompanied by a maximum at the zone center for other correlators. The oscillations 
indicate alternating orientational fluctuations, while the maxima at the zone center originate from 
nematic-like orientational fluctuations. For a < 2.5 and b < 0.35, the oscillations are weaker, leading 
to no marked maxima at the Brillouin zone edge. For a > 3.0 and b < 0.35, no oscillations occur any 
longer. For many of the orientational correlators in reciprocal space, an increase of a at fixed 6 or 
vice versa leads to a divergence at the zone center q = 0, consistent with the formation of nematic- 
like long range fluctuations, and for some oblate and prolate systems with b < 1.0 a simultaneous 
tendency to divergence of few other correlators at the zone edge is observed. Comparison of the 
orientational correlators with those from MC simulations shows satisfactory agreement. From these 
simulations we also obtain a phase boundary in the a — b-plane for order-disorder transitions. 

PACS numbers: 61.43.-j, 64.70.Kb 



I. INTRODUCTION 

The experimental, numerical and analytical study of 
structural properties of simple liquids is a well estab- 
lished discipline of condensed matter physics. In the 
center of such investigations is the static structure factor 
S{q). There are powerful integral equations allowing an 
approximate calculation of S{q) 1]. The starting point 
is the Ornstein-Zernike (OZ) equation, relating the to- 
tal correlation function h{q) and the direct correlation 
function c{q). An additional closure relation, like the 
Percus-Yevick (PY) approximation, then allows to de- 
termine h{q)^ from which S{q) follows from 



S{q) = l + ph{q), 



(1) 



where p is the number density of the liquid. Applica- 
tion of the PY approximation to a liquid of hard spheres 
yields good agreement with the exact result for interme- 
diate values of p 1]. However, the crystallization of hard 
spheres cannot be described by PY theory. 

The extension of the OZ equation and the PY approx- 
imation (or other closure relations) to molecular liquids 
is straightforward P, Q and has been applied extensively 
(see, e.g., 0,0). As for simple liquids, PY theory usually 
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does not yield an order-disorder phase transition. There- 
fore, it was quite surprising that a recent application of 
the molecular version of that theory to a liquid of hard 
ellipsoids of revolution with aspect ratio Xq has allowed 
the location of a phase boundary in the p — Xp-plane, at 
which a transition to a nematic phase takes place 

Much less analytical work exists for molecular crys- 
tals [Ml^. These are crystalline materials with, e.g., a 
molecule at each lattice site. One of the main interests 
concerns phase transitions of the translational and rota- 
tional degrees of freedom (see, e.g., the review JJ). These 
transitions are influenced by the translation-rotation cou- 
pling But one also expects phase transitions if the 
crystal is assumed to be rigid. For instance, let us fix 
hard ellipsoids of revolution with their centers at the 
lattice sites of, e.g., a simple cubic lattice with lattice 
constant equal to one. With a and h the lengths of the 
rotation axis and the perpendicular axes of the ellipsoids, 
respectively, an increase of a for fixed b may yield at a 
critical length ac{h) a transition to an orientationally or- 
dered phase, because the entropy of the ordered phase is 
larger than that of the disordered phase. 

Besides phase transitions, the study of the orienta- 
tional structure of molecules on a rigid lattice is of inter- 
est, too. Changing temperature will influence the steric 
hindrance between the molecules. The same happens 
for hard ellipsoids when changing a and h. Analogous 
to simple and molecular liquids, one can quantify such 
static orientational properties bythe one-particle orien- 
tational distribution {n) d E H ES 111 111 and by 
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the orientational correlation function where 
the orientation can be characterized, e.g., by the Eu- 
ler angles (0, 0, x) or, for axially symmetric particles, 
by {(f>,0). Then the following questions arise: How to 
compute p^^-* (f2) and, above all, the correlation function 
Gnn' by an analytical method? Docs the result 

for G„n'(ri, ri') allow to locate a phase boundary where 
a transition to an orientationally ordered phase occurs? 
To provide answers to these questions is the main moti- 
vation of our contribution. 

Our paper is organized as follows. In Sec. m we will 
introduce the model and the basic physical quantities like 
the one-particle distribution function p^^'^{^) and the ori- 
entational correlation function Gnn' (0,17'). The analyt- 
ical approach of calculating G„ri'(0, SI') or its transform 
Sxx' (q) from the OZ equation in combination with the 
PY approximation is described in Sec. IIIII Results from 
PY theory for hard ellipsoids of revolution on a simple 
cubic lattice will be shown in Sec. II VI and compared with 
those from MC simulations. This section also presents 
a short discussion of the phase transition for the ellip- 
soids from an orientationally disordered to an ordered 
phase. The final Sec. contains a discussion of the re- 
sults and some conclusions. We add a couple of appen- 
dices, including extensive technical manipulations needed 

in Secs. rrmiTvi 

II. DISTRIBUTION AND CORRELATION 
FUNCTIONS AND THEIR SYMMETRIES 

A. One- and two-particle distribution and 
correlation functions 

We consider a three-dimensional periodic lattice with 
A'' lattice sites and periodic boundary conditions. If the 
n-th lattice site has the position x„, the difference be- 
tween two sites is the vector x„„' = x„/ — x„. We as- 
sume the lattice to be rigid with lattice constant equal 
to one. At each lattice site we fix a rigid molecule, not 
necessarily with its center of mass. Restricting to lin- 
ear molecules, the orientation of the molecule at site n is 
given by r2„ — {(j)n,On) & S"^. The third Euler angle x„ is 
irrelevant for our purposes. Extension of our theoretical 
approach below to arbitrary molecules is straightforward. 
The interaction energy y({J7„}) is assumed to be pair- 
wise and the classical Hamiltonian is given by 

^ 1 

7J({0„}, {!„}) = Yl 2 1"!"' + ya^n}) , (2) 

n=l 

where 1„ and I(S1„), respectively, are the angular momen- 
tum and the tensor of inertia of the molecule at site n 
in the space fixed frame. Since we will investigate static 
quantities only, the kinetic part of i?({r2„}, {!„}) does 
not matter. 

In order to describe the orientational degrees of free- 
dom, we introduce the microscopic one particle density 



p„(ri), O €E S*^, at lattice site n and its associated fluctu- 
ation defined by 

Pnin) = sin\nn) , (3a) 
Spnin) = pnin) - (p„(r!)) , (3b) 

where 5(0|0') = sinO 6{9- e')5{<l)- (J)'), and ((•)) denotes 
canonical averaging with respect to if({r2„}, {1„}). Note 
that the f2„-dependence of p„ and 6pn is suppressed. The 
one-particle orientational distribution function is given 
by 

p«(0) = (p„(0)), (4) 

which is n-independent due to the lattice translational 
invariance of H{{Qn}, {In}), and the two-particle distri- 
bution is defined as 

p^„),(0,r!')EEp(f) ,(0,r!') = (p„(r!)p„,(r!')), n^n'. 

(5) 

(2) 

Note that p^^, is defined for n ^ n' , only, and that the 
index pair nn' is fully equivalent to the vector x„„/ . Mak- 
ing use of Eqs. it follows 

J p^^\fi)dn = 1, (6) 

y'pi^j,(0,0')dO = pW(0') (n^n'), (7a) 

J p5 {n, O') dO' = (0) (n ^ n') . (7b) 

It is important to realize that in contrast to isotropic 
molecular liquids the one-particle distribution p^^^ (O) de- 
pends on Q. An additional important feature may occur 
for particles with hard core interaction, e.g. hard ellip- 
soids, which are big enough. In that case there may exist 
a non-empty subset K of such that p^^^(f2) = for 
all n G K. Hence, the orientation of each particle is 
restricted to the complement K = \ K C S^. In 
the following, we need to distinguish between the two 
cases K — and K C S'^, where K is the area on 
which can be reached by particle orientations. Clearly, 

it is p^^l{n,n') = 0, if p^^Hn) = o or p^^Hn') = o, but 

not vice versa. Three examples for p*-^) (0) are shown in 
Fig. n obtained from MC simulations for hard ellipsoids 
on a simple cubic lattice. 

Next, we introduce the orientational density-density 
correlation function G„„'(S1, fl'). It is the correlation of 
the fluctuations of Pn{^) at lattice sites n and n': 

Gnwi^M) = {5pnm5p,,,{n')) . (8) 

By use of Eqs. ®-© we get 

G„„. (o, n') ^ ,5„„,p(i) (0) 5[n\n') - (0)^(1) (o') 

+ (l-<5„„0p5(f^,^^')- (9) 
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FIG. 1: MC results for p'-^^il) for hard ellipsoids of revolution with a = 3.6, 6 = 0.24 (left), a = 1.2, b = 0.88 (middle) and 
a = 0.56, 6 = 1.1 (right) on a sc lattice. Orientations on obtained from MC runs are projected along the fourfold lattice 
direction. Parts of K occur along the two-, three- and fourfold lattice directions, depending approximately on whether 2J± 
exceeds the site-site spacing along one of these directions. Some circles approximating the edges of the parts of K are shown 
as an aid for the eye. 



Gnn'i^, = Ol^lin, n') + Gl^l{n, n') consists of a self 
part and a distinct part, which are explicitely 



(lOa) 

Gi'l{n,n') = (1 - <5„„,) {p^^l{n,n') - p('H^)p^'\^')) 

(lOb) 

Due to the properties ^ and of the particle distri- 
bution functions, Gnn'i^, fulfiUs for aU im' 

J Gnn'{n,n')dn^ j G^n'{n,n')dn' = q. (ii) 



S2 



S2 



For the lattice system, the pair and total correlation 
functions (7„„' (17,17') and hnn'i^^fjD are introduced in 
the same manner as for a liquid [l|, Q] : 



{n ^ n') , (12) 



p(i)(17)p(i)(f]') 

hnn'{n,n')=gnn'{flM')-'^ {n ^ u') . (13) 



Here, some comments are in order. First, g„„' and /i„„' 
are defined for n ^ n' only, and if /sT C 5^, for 17, 17' G if 
only. Second, /i„„'(17, 17') in general obeys 



;i„„.(17,17')rfl7 7^ 0, y /i„„-(17,17')dl7' 7^ {n^n') 

(14) 

The same is true for (7„„'(17, 17'). Third, due to 
Eqs. ©, O, linji and (HSJ it is 



p(^)(17)/i„„'(17,17')(il7 = {n^n'), (15a) 



S2 



/i„„/(17,17')p(i)(17')dl7' = (riT^n'). (15b) 



S2 



Last, in the asymptotic limit of large particle separations, 
g„„'(17, 17') and /i„„'(f2, 17') behave like 



lim 5„„'(17,17') = 1 , 

lim hnn' (^7 ^0 = ! 



(16) 
(17) 
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independent of the direction of x„„/ . This follows due to 
lim|,,^,^,|_o, ^ in Ml agree- 

ment with the behaviour of a liquid system. 

Making use of Eqs. ifTUI) . l(T^ . ifT^ and introducing 



D{n, n') = 47r (5(17|J7') - (n) (17')) 

we can rewrite Gl"J,{fl, Q'), a = s, d, as follows: 



(18) 



G[^l,{n,n') ^^Snn'D{n,n'), (i9a) 

(19b) 

As for molecular hquids 0] (see also, e.g., 
it is useful to expand all orientation-dependent functions 
with respect to a complete set of functions, determined 
by the rotational symmetry. If V is the point symmetry 
group of the lattice and Vm the symmetry group of the 
molecules, one can use basis functions for irreducible rep- 
resentations of the symmetry group of p(^^(f7), which is 
a subgroup of V ® Vm 0, M, EOI ■ For axially symmetric 
particles, these are linear combinations of the spherical 
harmonics Yx{n), A (Im) HQ. On the other hand, 
the spherical harmonics itself can be taken. We have 
chosen the latter possibility in order to keep similarity to 
molecular liquids. Consequently, we have for any func- 
tions /(ri) and F„„/(r2, ri') their A-transforms and the 
corresponding inverse transformations: 



(20a) 
(20b) 



h = j f{n)Y^{n)dn, 
fm = ^(-z)-'/,n(f7). 



A 



(21a) 



AA' 



(21b) 



The purely imaginary prefactors in H20a|l and (|21a|l are 
taken for technical convenience. 

Finally, we can use the lattice Fourier transform due 
to the lattice translational invariance. It is restricted to 
the first Brillouin zone of volume Vbz- For example, the 
transform of the site-site matrix elements (|21af) and its 
inverse are given 

^^AA'(q) = ^e''i'''".'i^w,AA', (22a) 
Fnn',xx' = / i^AA'(q)e~"'"""'d'<7. (22b) 

l.BZ 



B. Symmetry relations for the one-particle 
distribution and the correlation functions 

Let us first discuss the one-particle orientational dis- 
tribution function p'''^\^l). p'^^\^l) must carry the full 
point symmetry V of the underlying periodic lattice, and 
also the symmetry Vm of the particles Neglecting 
the latter for a moment, p'-^-'(f2) for axially symmetric 
particles can be expanded into a series of all T'-invariant 
combinations Y/n, (f2) of spherical harmonics. This ex- 
pansion reads (cf. pObfl l 



Ini 



(23) 



Here, the numbers n; are the multiplicities of the unity 
irreducible representation contained in the point group 
representation of V established by all spherical harmonics 
of order I. Then, by the invariance requirement under 
symmetry operations of the particles, H23|) can eventually 
be further simplified d, . 

If the lattice is cubic, V = O/i, and for any kind of 
axially symmetric particles we have 



,(1) 



+p^i^%i{n) + o{i = iQ). 



(24) 



The Oh cubic invariants Ym{VL) = yoo(f^) = (47r)-i/2, 
Yji jO), Yfii (r^) and 181(0) are real functions and given 
in [ill, lij up to factors {At:)-'^/'^ p~K 

In Sec. IIVI the canonical averages (1a) = 
P^^K^) >a(0) dn are needed. The values (Kn), {%{) 
and (Fgi) have been calculated by MC simulations (see 
Sec. lIV|) for several values of a and 6, and the nonvanish- 
ing (Y\) up to Z = 9 are 



(i^oo) = 


1 

V47r' 


{y4o) 




(n±4) = 


730 . 

12 ^^^'>' 


(>60) 




(i^6±4> - 




(Yso) 


- f (ft,): 


{Ys±i) = 


^(Y ) 
24 




- 48 <*^i" 



(25) 



Next, we investigate the symmetries of G„„' (O, il'). 
From the definition ((HJ it follows immediately 



Gnn'(0,57') — G'„'n(0', O). 



(26a) 



If the inversion / belongs to V, the inverted Gnn'{^, O') 
must match the old one by use of Eq. ijHJl, since 
H{{fln}, {In}) remains unchanged under inversion then: 



G„'„( — O, — O') — Gnn' (O, O'), 



(26b) 
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where = {9, 0) < — > — f2 = 7J7 = (tt — 0, (/j + tt) has been 
used. If the symmetry group V contains rotations R, the 
rotated Gnn'i^,^') must be the original one: 



(26c) 



The knowlegde of the transforms on the Ihs of Ij26b|l 
and H26c|l is sufficient to determine the transform 
Ggx„„, (ff^j 5^') for an arbitrary point symmetry oper- 
ation g, which must be Gx^^, (O, ^i') if g ^ V. 

The properties H26(l affect the matrix elements Gnn',\\' 
the following way: 



Gnn'.XX' 
Gnn'.XX' 



^ n' n,X' X 
(-1) 



Gr, 



.XX' 



n"ni"' 

X G 



m' tn 



,„{R) 



(27a) 
(27b) 

(27c) 



In H27c|l . Wigner's generalized spherical harmonics (rota- 
tion matrices) are used 01 • To calculate the correct rota- 
tion matrix elements for (|27c|l . one uses the three Euler 
angles carrying some coordinate frame, in which x„„/ is 
assumed to be fixed, into a new, symmetry-equivalent one 
such that the rotated vector x„„/ coincides with i?x„„/. 
The behaviour of {Y/m(17)} under complex conjugation 
yields a fourth property: 



Gn 



Am A' 



(-1) 



.l — m.l' — r 



(27d) 



Eqs. H27(l are translated to the Fourier transformed matri- 
ces S(q) = 47rG(q) = 47r(GAA'(q)) (see Eqs. (EZj), (EHJ) 
due to 



5'AA'(q) 



(28a) 

^AA'(q) = (-1)'+' 5AA'(-q), (28b) 



X S';m",i'm"'(q) , (28c) 
SlmMra'{ci) = (- 1)'+''+"+"' 5r_„,,,_,„, (-q) (28d) 



?a|l demonstrates that S(q) is Hermitian, but G„„' in 
general is not. 

The symmetry of the particles can bring about 
extra characteristics of the matrix elements |^ 0. 



For axially symmetric particles with inversion symme- 
try, Git\(n,fi') = Gill (-n, 17') = Gitl (17, -17') = 
^in'(~^'~^') vaHd, and the self part fulfills 
G''nl{n,VL') = Gl"2(-f7,-f7'). Consequently, Gi"2 can 
have nonzero elements for I and V even or I and V odd, 
respectively, while G^^'^', has nontrivial elements for I and 
V even, only. 

We want to conclude this section with the remark for 
cubic lattices, that by symmetry the knowledge of the 
correlations for only ^ of all lattice vectors or ^ of the 
volume of the first Brillouin zone is necessary to calculate 
the correlations for the complete lattice or Brillouin zone. 

III. ORNSTEIN-ZERNIKE EQUATION AND 
FERGUS- YEVICK APPROXIMATION 



Similarly to simple and molecular liquids, we will in- 
troduce the direct correlation function Cnn'(J7jf7'), which 
is related to /i„„/(ri, fi') by the OZ equation Since 
c„„'(il, f2') is determined by the inverse functions of 

^m!'(^'^') G'„„'(ri, f2'), one has to be careful be- 
cause of relations which imply that a constant func- 
tion /(ri) = const is an eigenfunction of g[j2'(^^j and 
G„„'(ri, ri') with eigenvalue zero. Therefore, these in- 
verses do not exist on the one-dimensional subspace of 
constant functions. This is a new feature occuring for 
rigid molecular crystals. The solution of this problem 
is easy using the A-transforms -^Snn' D\\i and Gnn' ,xx' 

of g\IIi{^t^') and G„„' (Jl, i7'), respectively. Because of 
Eqs. it follows that the first rows and columns of 

the matrices D and G„„' are zero, i.e. 



Dxx' = , G 



7l7l' ,X\' 



= 0, ; = and/or = 0. (29) 



This is reasonable since quantities with / — and/or 
?' = do not describe orientational degrees of free- 
dom and therefore are unphysical. Obviously inversion 
must be restricted to the A-transform with II' > 0, 
which is equivalent to restrict inversion of G^^*,', (fJ, fi') 
and G„„'(r2, fi') to the space of functions having no con- 
stant part. Since an arbitrary function (fJ, f2'), and 
maybe c„„/(0,r2') as well, need not fulfill relations ltTT|l . 
we associate with any function (fi, 17') on 5*^ x S'^ 
the function 



F°„, (17, 17') = Fnn' (17, 17') - ^ y Fr,n' (17, 17') dn-^J F„„, (17, 17') dl7' + J J F^^, (17, 17') dfl dl7', (30a 



which can be rewritten as 



F°„,(17,17') ^ JJ (^^(17|17") - i-^ F„„-(17",17"') (^,5(17"' |17') - ^ ) dVl'dVl"'. (30b) 
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It is important to note that the two functions Fnn' 
and F°^,{Q,{1') differ at most by the sum of a function 
of n only, a function of fi' only and a constant func- 
tion, and that (17, f]') fulfills relations l(TT|l . In gen- 
eral, Fnn' cannot be recalculated from F°^, (fi, fi'), 
since l|Snjl describes a projection. The associated projec- 
tion operator R = = is 



R{n, n') = 5(n\n') = 5°(n\n') . (sia) 

47r 

By Eq. (jSObp . we can present F°^,{il,n') in the short- 
hand notation 



F:n'{^,^') = iRPnn'Rm,n') 



(32a) 



The elimination of the first column and row (cf. Eq. (|29|) ') 
of the A-transform Fnn' .w can be done similarly by use 
of the A-transform 

Rw = j''"' (^AA' - '5a,oo<5oo,a') (Elb) 
of i?(f7,ri'). This leads to 

Pnn',\\' — ^ R\\"Fnn',\"\"'R\"'\' (|22b) 
A" A'" 

Relations imply that 

G„„,(f7,l]') = G°„,(f7,f7'), /?(r!,r!') = i?°(f7,f7') (33) 
for all nn' , but, in general, 

hln'{n,n')^hnn'{n,n') {n^n'). (34) 



For later purposes it will be helpful to rewrite Eq. Ijl9bp 
replacing /i„„/(r2, fi') by (fi, fi'). Eqs. l(T^ can be 
used to show that the rhs of Eq. Ijl9b|) is indeed of the 
form o and that the function D(r2, fi') can be introduced 
as follows: 

p'^^\n)hnn'{^M)P^^\^') 
1 



(47r)2 



D{^, n") hnn'{^\ ^"') D{n"', n') dn" dn'" 



(47r)2 



D{n, ii") Kn'{^\ ^"') D{n"', n') dn" dn'" 



{n ^ n) , 



(35) 



where the last line involves {n, n') instead 

of hnn'{n,n') since D{n,n') = {RD){n,n') = 

{DR){n,n') = {RDR){n,n') projects out the constant 
parts of hnn' i^, ^')- Contrary to the general case, it will 
be shown in Add. 1X1 that h'^,^,{n, n') already determines 
hnn'{n,n') uniquely (cf. Eq. HA4p '). The last step is just 
to write down G",^, (il, n') in the new form 

G°„,(f7,17') = j-6nn'D{n,n') + il-6nn') 

X j j D{n,n")h°nn'{.^" M"')D{n"\n')dn" dn"' (36) 



The tensorial static structure factors S\\i{q) = 
S'^y{q) also have their first columns {V = 0) and rows 
{I = 0) vanishing and can be defined by the two equiva- 
lent expressions 

5aA' (q) = ^ {5pl (q) Spy (q)) = 4^ Gaa' (q) , (37) 

where the factor 47r is introduced as the inverse of the di- 
rection average of the one-particle density p'^'^\n). The 
transform of Eq. (|36(l leads directly to the concise ex- 
pression for the structure factors Sw (q) in matrix form. 



S(q)=D + i-Dh°(q)D, 



(38) 



where for the Fourier transform h° (q) it is assumed that 
hnn{n,n') = 0. Eq. closely resembles the corre- 
sponding relations for simple Ql and molecular liquids 0] ■ 
There is a second problem occuring if hard core inter- 
actions lead to K C S^. As already discussed in Sec. Ill Al 
p^^^(ri) and therefore G„„'(ri, Jl') vanish in that case on 
the subset K C for all nn' . Therefore, any function 
g{n) which vanishes on K but is nonzero on the com- 
plement K is an eigenfunction of G„„'(ri, n') with eigen- 
value zero, too. This problem can be properly solved as 
well. It is reasonable to discuss hard core interactions 
leading to if C S*^ separately. Therefore, the next sub- 
section nil Al will present the OZ equation and the PY 
approximation for K — S^, and in subsection IIII Bl sys- 
tems with K C will be treated. 



A. Molecular crystals characterized by K — 

Of course, K — for soft potentials. But also for 
hard core interactions it can he K — S^, provided size 
and shape of the particles are properly chosen. This can 
be the case, e.g., for slightly aspherical hard particles. 

Above, everything has been reduced to functions of the 
form o, and also the direct correlation function can only 
be defined as c°„, (fi, il'), at least at this first stage of 
definition: 

cl,An,n') = 
(g(^)°) {n,n') - ({G°r') {n,n') (39) 

The reader should note that, in contrast to pl^J^,{n,n') 
and h'^^,{n,n'), the self part of the direct correlation 
function exists and is well defined by Eq. (|39|l . 

The function G°„, (ri,il') has an inverse on the sub- 
space of functions g°{n) only, since Gnn'{n,n') = 
G°^^,{n,n') = {RGnn'R){n,n') maps constant functions 
to zero on S^. This inverse is determined itself only in 
the form o and fulfills 



E 



52 



(G°)-^)° ^{n,n")G°nnn'{^'\^')dn" 



6nn'R{n,n'), n,n'es^ (40) 
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Note that the rhs of Eq. ^ does not involve S{n\n') but R{n\n') = S°in\n'). 



An analogous equation defines the inverse of the self part Eq. HlUa|) . which is, contrary to the full inverse above, 
known explicitely: 

\ 52 / 

The derivation of the lattice OZ equation causes no problems having the tools presented in Sees. ITU and IIIII above. 
The procedure is to substitute ((G'°)-i)°„, (^2, rj') from Eq. ^ and G°„,(r2,rj') from Eq. §^ into Eq. gUJ, at the 
same time using that G^'^, 17') = ^Snn'D{n,n') and fg^{{D°)-^)° {n,n") D{n" ,n') dfl" = R{n,n'). Then, the 
Ihs of Eq. becomes 



^^„..„.D(f}", 12') + (1 - <5„,,„o y J D{n", ^"') K„^,{n"', n"") D{n"", n') dn'" dn"" | dn" 
= s^n'R{n,n')-^ J C^,in,n")D{n",n')dn" + ^{i~s^n') / hi„,{n,n")D{n",n')dn" 
-^^(i-<5„»„o c°„4^7,^7'')i?(^7'',^7''')C'„'(^^'^^^'''0^(^^''^^^0^^^^''^^^^^^ (42) 

for 12, 17' e S*^, which must match the rhs of Eq. H4()|l . i.e. 5„„/i?(12, 12'). (5„„'i?(12, 12') cancels. Operating with 
Jg2 ■ ■ ■ ((-D°)^^)°(r2', 12""') dr2' on both sides gives the OZ equation for the lattice in the site-site angular representation 
for 17,17' e S^: 

{i-s,,r.')K„A^M)^ci^'{^M) + ^ E j j cl,A^,n")D{n" M")hl„A^"' M)dn" dn'" . (43a) 



The associated A- and Fourier transformed matrix equa- 
tion, which has the most significant form and which is 
used for the numerical work, reads 

h°(q)-c°(q) + -^c°(q)Dh°(q), 62b) 

47r 

where again it is assumed that /i„„(f2, 12') = 0. 

This result is almost the same as for molecular liq- 
uids 0- There are two main differences. First, we have 
to use all matrices in their form o, i.e. the first row and 
column of the original matrices are skipped, since they 
are zero. Second, there appears the matrix D, which is 
identical to G^*-', apart from a prefactor, and third the 
Fourier backtransform of h° (q) has to fulfill ft.„„ aa' = 
for all AA'. 

Up to now, the equations are not closed. Since nei- 
ther the total correlation nor the direct correlation func- 
tion is given, the previous concepts are almost useless if 
one is interested in an analytical approach to determine 
the structure factors H37|) . (|38|l . An additional equation, 
called the closure relation, must be used to find a self- 



consistent solution for /i„„'(r2, 17') and c„„'(17, 17'), as for 
simple and molecular liquids. It has been our intention 
to follow as close as possible the established lines of liquid 
theory, and so we chose the most straightforward analo- 
gon of the PY approximation for the lattice, which is for 
n ^ n' 

c„„, (17, 17') = /„„, (17, 17') (.9„„, (17, 17') - c„„. (17, 17')) 

= fnn' 17') (1 + Ku' 17') - Cnn' i^, f^')) ■ 

(44) 

Note that Eq. ^ involves c°„, (17, 17'), only, but in ^ 
the full function c„„'(17, 17') appears. In Add. IdI it will 
be shown for hard particles and n ^ n' that c°„, (17,17') 
determines c„„' (17, 17') uniquely. /„„' (17, 17') is the Mayer 
/-function, which is for n ^ n' 

/„„,(17,17') = exp{-/3K„'(f^n,^^«')} - 1- (45) 

For hard particles, the pair potential is Vnn' (1^, = 
(if the pair (nl7,n'17') has no overlap) or T^„„' (17, 17') = 



8 



oo (if the pair has overlap). This imphes that all static 
properties are athermal and that it is 



, no overlap , 
— 1 , overlap . 



(46) 



The range of the /-function for hard ellipsoids is 
max(a, 6), if the ellipsoids are fixed with their centers 



of mass on the lattice. In accordance with the theory 
of liquids of hard particles, lEjl yields (/„„' (fi, fi') = 

(2) 

< — > (n, ri') = while c„„'(ri, f2') remains undeter- 
mined, if (nfi, n'fl') has overlap and c„„' (SI, fl') = while 
f;„„/(f2, i7') remains undetermined, if {nil,,n'D,') has no 
overlap [IQ. 



I 

Eq. (|44|l should be available in its A-transform, since the numerical solution of the OZ equation using the PY 
approximation is most conveniently done in terms of matrix elements. From Eq. H44|l . it is straightforward to deduce 
that 



l+l" -V" +1"" -I'"" 



{21" + 1){21 



III! 



1) 



47r {21 + 1) 



(2r + l)(2r"" + l) 



At:{2V + 1) 



Cnn', AA' 

A" A'" A"" A" 

C{l"ri, 000) C{l"'l""'l', 000) C{l"ri, m"m""m) C{l'l"'l"'\ m"'m""'m') fnn'.yy {9n 
where the Clebsch Gordan coefficients C {111213,1^11111211113) a re given in f^- 



.A""A'' 



^7in',A""A'"" ) ; 

(47) 



The numerical use of Eq. 147() requires the knowledge of 
the matrix elements /„„',aa' and also5„„'^AA' or hnn'.xx', 
n ^ v! . These are derived in Apps. IbI and lO Having 
calculated the matrix c„„' then, the Fourier transform 
of c°„, = Rc„„'R serves as input for the OZ equa- 
tion iPEj) . 

The fact that the OZ equation relates c°„,(ri, J7') and 
/i°„,(17,fJ') (for n ^ n') to each other but the PY ap- 
proximation involves c„„'(n, J7') and /i„„'(f2, il') requires 
some discussion. In App. ^we prove that /i„„/ (fi, $7') 
is uniquely determined by ($7, ^2') (cf. Eq. ^^). 
Therefore, having determined /i"^^/ (O, J7') from the OZ 
equation, Eq. (|A4(I yields /i„„'(f2,n') (n ^ n') and this 
in turn yields c°„, (^1,17') by use of the PY approxima- 
tion . In order to calculate Gnn' ,\\' or G'AA'(q) one 
also needs the self part of the direct correlation function. 
Taking n = n' in Eq. H43a(l yields the matrix equation 



1 

47r 



Dh° 



(48) 



i.e. the self part c°„(r2, 17') is determined by the distinct 
parts of c°„,(n,f]') and /i°„,(f2, J]')> only. This feature 
exhibits a further difference to the OZ equation and PY 
approximation for liquids. If if C S"^ , additional new fea- 
tures occur which will be discussed in the next subsection 
and App. [E] 



B. Molecular crystals characterized by K C 

As already stressed in Sec. Ill Al there may exist a non- 
empty subset K of on which p'^^^(O) and Gnn'{^, 
vanish. This case only occurs for hard particles, if their 
size exceeds a certain limit. Accordingly, any function 
/(ri) or F„„/(r2, fi') must be restricted to e K, 



denoted by f^{Q) and F^^^,{Q,n'), respectively. It is 
easy to prove that Eqs. (|ll|l read now for all nn': 

Gl, {n, n') dn= f g^, {q, n') dn' = o. (49) 



K K 

This suggests to introduce for each function i^^^, (fi, U!) 

f5°(17,17') = {R^F;^^,R^){^IM) (50) 
with the projector 



, otherwise 

having the A-transform 

Y^{n)Yy{n) dn 



(51a) 



K 



1 



^— / Y;:{n)dn / Yy{n')dn', (5ib) 

K K 

which are generalizations of Eqs. (|31|l and ^A2\ . Note that 
I if I is the total space angle covered by the area K C S^, 
and that in Eqs. Yx{n) and Yx'{n') have to 

be restricted to K. li K C S'^, {Yim{^)} is no longer 
an orthonormal set. The latter implies the very impor- 
tant fact that if if C 5^ there exists no angular function 
F^{n,W) for which the associated matrix F^, is the 
unity matrix. Therefore, inversion ofF^,;?{n,n') and of 
^nn' iTiust be performed with respect to {il, fl') and 
R-^, respectively: 



E 



'F 



K.i 



K,< 



{n,fi")F^;°^,{n",n')dn" 



K 



5nn'R'Hn,n') 



(52) 
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and similar for the A-transforms. Taking these modifi- 
cations into account one can follow each step done in 
subsection IIII Al Finally one obtains the OZ equation 
for the case K C which follows from Eq. H43a|) by re- 
placement of S'^ and o by K and K, o, respectively. The 
factor however, remains. The same replacement has 
to be done for all functions occuring in the PY approxi- 
mation igni: (03 • 

So far the OZ equation and the PY approximation for 
K C are directly related to the corresponding equa- 
tions for X = S*^ by the replacements mentioned above. 
A new feature is related to solving the OZ equation 



(q) 



(q) 



47r 



c^^^°(q)Dh^-°(q), (53) 



for h^'°(q), which is needed for the self-consistent solu- 
tion. In the case K — S^, one performes a simple matrix 
inversion. For K C S^, this can be done in the following 
way. First, rewrite Eq. (|53|l as 



R''-^c^'°(q)D)h^-°(q)=c 



1 



K,o 



(q). (54) 



Note that we have to use R^h^^'°(q) = h^'°(q), since 
there exists no function F^{Q, fi') with F^, = Sw . The 
first rows and columns of the matrices and — 
|^c^'°(q)D vanish. The inverse of the latter matrix 
with respect to R-'^' is given by R^ + -^h^'°(q) D, due 
to Eq. Then it follows 



h^'°(q) = R^ c^'°(q)D 



An 



c^'°(q): 



(55a) 



where the lower index R^ indicates inversion with re- 
spect to R^. In the next lines, 1 means the full unity 
matrix or the unity matrix for I, I' > 0. It is interest- 
ing that the matrix 1 — |^c-^'°(q)D still has an in- 
verse with respect to the unity matrix 1. This inverse 
is 1 -I- |^h^'°(q) D, according to Eq. H53|) . Consequently, 
an equivalent form to H55a|l is 



q) 



-K.o 



47r 



(q)Dj^ c 



(q) , (55b) 



Indices K or o for the inverse on the rhs of Eq. (|55b|l have 
been left out, since it describes no angular function like 
the inverse in (|55a|l . It is easy to prove that the matrices 
h^'°(q) calculated due to Eqs. (|55|l are hermitian, if D 
and c^'°(q) are. 

Because of these complications for K C S^, on one side 
one must solve the OZ equation using PY theory on the 
restricted angular space, i.e ft, G K. But on the other 
side all physical quantities entering the OZ and PY equa- 
tions as input are continuous under variation of the pair 
potential y„„' (fi, $7'). Even though the physical part K 
of may experience a jump from K = to K C S^, 



the physically relevant structure factors H37|) , (|38(l should 
behave continuously, too, if we work on S^, assuming fi- 
nite interaction V < oo whenever in reality we have hard 
core overlap, and then turn on the hard cores by tak- 
ing the limit V — > oo. In App. ^ it will be shown 
that if hard core interactions lead to if C S*^, one can 
indeed determine solutions h'^^,{n,il') and c^, (il,f2') 
of the OZ/PY equations for n,n' e and V = oo, 
which lead to the same structure factors as the functions 
/i^„,(ri, r?') and c^, (fi, r2'), introduced above. Conse- 
quently, for practical purposes it is much easier to solve 
the OZ and PY equations defined for fi, W G in 
Sec. IIII even for if C S"^ to calculate the strucure fac- 
tors S'AA'(q)- 



IV. RESULTS FOR HARD ELLIPSOIDS ON A 
SC-LATTICE 

The self consistent solution of the OZ equation and the 
PY approximation has been done numerically. In order 
to check the quality of these solutions we have performed 
MC simulations, which also allow to determine the phase 
boundary between orientationally ordered and disordered 
phases. We stress that the investigation of the phase 
transition has not been our major motivation. Therefore 
we have not attempted to verify the phase boundary by 
more complicated MC algorithms. Before we describe the 
numerical solution of the OZ/PY equations in Sec. lIVBl 
let us present some details of the MC simulations in the 
first subsection. The results from both approaches will 
be discussed in Sec. IIVCI 

For both the numerical solution of the OZ equation us- 
ing the PY approximation and the MC studies the over- 
lap criterion of Vieillard-Baron for hard ellipsoids of 
revolution is used. 

The numerical studies are done in terms of matrix el- 
ements for all A up to a certain limit Zmax, refering to all 
correlation function matrices up to the /-function ma- 
trices, for which the maximum Z-value is 2 /max- Since 
the ellipsoids have inversion symmetry and are assumed 
to be fixed with their centers of mass on the lattice, the 
matrices G„„', h„„', c„„/ and f„n'j n ^ n', have non- 
vanishing entries for I and I' even only, whereas the ma- 
trix D, and therefore also G„„, has also non- vanishing 
entries for I and I' odd. Despite of these {W odd)-entries, 
the OZ equation H43b|) reduces to an equation for the 
blocks LI' > 2 and even. 



A. MC simulations 

The MC simulations are performed for systems of 
16 X 16 X 16 simple cubic lattice sites and periodic bound- 
ary conditions. For systems having very long-ranged cor- 
relations along a certain lattice direction, test runs with 
32 X 32 X 32 particles are done, but the results are only 
slightly different and not shown in this work. 
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The only possible move is a rotation for a randomly 
chosen particle perpendicular to its orientation axis u„ 
by an angle < < 0,nax < tt, where cos 9 is at random, 
and a subsequent rotation with respect to its original ori- 
entation by a random angle < (/) < 27r. If this move 
leads not to an overlap of ellipsoids, it is accepted, other- 
wise rejected, in which case the next move is not tried for 
the same particle, but for a new randomly chosen one. 

As starting configuration for each MC run parallel el- 
lipsoids are chosen. Each particle is moved on an average 
of 1000 times with ^^ax = f to get a rapid convergence 
to the disordered phase having a cubic symmetry, if pos- 
sible. This process can be followed measuring the biggest 
eigenvalues A+ of the Saupe tensor for the instanta- 
neous configurations at certain time values and also for 
subsystems of appropriate size. In the disordered, cubic 
phase the time average of A+ must decrease with the sys- 
tem size, since A+ is connected to the canonical averages 
(^2m)j which must vanish due to Sec. Ill Bl Next, each 
particle is moved on an average of further 1000 times, 
after each 10 times adjusting 6'max by ±0.02 to get a ac- 
ceptance rate of « 25%. Then, the system is equilibrated 
doing an average of 5 x 10* moves per particle, which is 
equivalent to about 2 x 10^ MC steps. After this, an 
average of further 10^ moves is done for each particle, af- 
ter each 10 moves storing A-|_ to calculate the correlation 
time r of 



A(t) 



(A+(t)A+)-(A+)(A+) 
(A+A+)-(A+)(A+) 



(56) 



In the production run, every t x 10 moves per par- 
ticle the current configuration is used to measure the 
values = i'HY*) (/ = 2, 4, 6, 8; to = 0, . . . , and 

pSUv ^ Pnl',xy ^ PxL,AA' (0 < x,y,z < 16; / = 
2,4; /' > < m < /; -m < to' < TO if ? = I'; -I' < 
to' < /' if Z 7^ /'). See Eqs. (PUajl . |[?Tal) for the definitions 
of p^^^ and p^xyzW- ^^^^ condition for the simula- 
tions, a standard deviation of lO^'' or less for the values 

P^xyz w been chosen. 

After the production run, a more serious evidence for 
the cubic symmetry than the behaviour of A+ is to check 
whether the measured data (Y\), have about the pro- 
portions given in Eqs. H25() and if the other values really 
"vanish". Then, the matrix d (see Eq. (IB2bll ) is calcu- 
lated by Clebsch-Gordan coupling '2*1 from the measured 
data and averaged over all cubic symmetry operations. 
After this, d is composed by the best estimates for the 
values (Kii), (^6i) and (Fsi), and its first row is used to 
extract these values. In the last step of evaluation, the 

(2) 

above missing matrix elements p\.y^ are calculated by 
symmetry. These are also subjected to all cubic symme- 
try operations, each referring a to point symmetry oper- 
ation on the complete system, and averaged. The third 
test for the cubic symmetry is then to check if these av- 
eraged correlations match the old ones. 

The resulting phase diagrams for prolate and oblate el- 
lipsoids are shown in Fig. |21 The thin solid lines charac- 



terize the closest packing of parallel ellipsoids. They rep- 
resent upper bounds for the phase boundaries for transi- 
tions to ordered phases with aligned ellipsoids. Whether 
there exist more complex, ordered phases, commensurate 
or incommensurate ones with even larger volume fraction 
than on the thin solid lines is not known. An interest- 
ing feature of these lines can be observed. For prolate 
and oblate ellipsoids there are characteristic pairs (a, b) 
at which cusps occur, indicating a maximum volume frac- 
tion. The light grey areas represent the transition region 
from an orientationally disordered to an ordered phase. 
The latter is not necessarily a phase of aligned ellipsoids. 
Transitions have been observed from the ordered to the 
disordered phases and also vice versa for systems with 
ellipsoids of small enough maximum linear dimension, 
having only few interaction partners. The hysteresis is 
small, indicating either a continuous or a weakly first or- 
der phase transition. The medium grey areas refer to the 
OZ/PY solution and are explained in the next Sec. lIVBl 
while the dark grey areas indicate overlap of the light 
and medium grey ones. 



B. Numerical solution of the OZ equation using 
the PY approximation 

In App. El we show that the correct static structure 
factors S\\'{q), even for K C S^, can be obtained by 
solving the OZ/PY equations for 17,17' S S'^. We have 
chosen this option since it avoids the use of the projectors 
and (cf. Sec lmBl and App.lQ. We also solved 
both equations for 17, 17' S iiT as described in Sec. IIII Bl 
but it turned out that the results are worse. The reason 
is probably the fact that the cut-off Imax, at which the 
matrix equations are truncated, leads to a rather crude 
approximation for and . For example, some of the 
truncated projectors are significantly non-idempotent. 

The numerical solution of the OZ/PY equations is per- 
formed by the iterative procedure described below usu- 
ally for lattices of size 32 x 32 x 32, periodic boundary 
conditions and /max = 4. The correlation function ma- 



trices Gxyz — Gy 



have, by symmetry and 



periodicity, only to be tabulated for 0<x<y<z< 16, 
and the same for the discrete Fourier transformed matri- 
ces S(q). 

For some systems characterized by (a, 6), we have 
solved the OZ/PY equations additionally for /max = 2 
and /max — 6. All the correlators we have investigated 
explicitely remain qualitatively unchanged, while the cor- 
relation lengths of these correlators become larger for in- 
creasing Imax- For some of the systems showing no con- 
vergence of the iteration scheme because of diverging cor- 
relation lengths for /max = 4 and 32 x 32 x 32 lattice sites, 
system sizes up to 128 x 128 x 128 have been used, but 
for none of these bigger system sizes convergence could 
be achieved. 

After each numerical inversion or Fourier transform the 
symmetries are checked to be within certain error mar- 
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a) prolate elliposids 




b) oblate ellipsoids 



0.7- 



0.4- 



0.3- 



0.2- 



0.1- 




1 1.1 1.2 1.3 1.4 1.5 1.6 



FIG. 2: Phase diagrams for a) prolate and b) oblate hard 
ellipsoids of revolution on a simple cubic lattice. Solid lines 
refer to the closest packing for parallel ellipsoids. Within 
the light grey areas an order-disorder phase transition occurs, 
while the medium grey areas indicate where the numerical 
solution of the OZ/PY equations starts to diverge. For a 
more detailed discussion see text of Sees . II V XI and HV Bl 



gins, and then each matrix is subjected to every symme- 
try operation under which it must be invariant. From all 
these symmetry-transformed matrices the average is cal- 
culated to rule out rounding errors to become too large. 

Before starting the iterative procedure, the input quan- 
tities are calculated: 

(a) For each pair (a, b) of ellipsoid axes, the /-function 
matrix elements //^/,„ in the r- frame |j| according to 
Eq. I|C5|I have to be calculated for all even I < I' < 2 Z^ax 
and < m < ma.x{l,l'). By Eqs. (|C1|) . the remaining el- 
ements are obtained. Then, these are transformed back 
to the laboratory frame. For the numerical integration a 



400 X 200 X 200-grid has been used. ^ 

(b) The values (I41), (^ei) and (Fsi) are taken from 
the MC simulation for the same pair (a, b) or an appropri- 
ate extrapolation beyond the MC cubic phase boundary, 
where needed. Then, Eqs. H25() yield the canonical av- 
erages of the spherical harmonics, which are coupled to 
give the matrices d and D, due to Eqs. (jB2bp and ljB2c|l . 
The iterative procedure is now as follows: 
(a) As input for the OZ equation in the n-th iteration 

step, we use the matrices ci^^ , which will be specified 
below. According to Eqs. (|31b|l and H32b|l . the matrices 



are obtained by just cancelling the first columns 



o,(n)' 



^xyz 

and rows. The matrix elements c^y^ 
lated for < x,y, z < 32 and subsequently Fourier trans- 



are then calcu- 



formed to yield Cp'^"\x' 



- c^^^' - 



■xyz 



. As very first input we choose 
otherwise. 



o,(n)' 



(b) The matrices 
calculated from 



in the n-th iteration step are 



pqr 



^o,(")' 



n > 1 . 



They generally will not fulfill the constraint h^^ 



(57) 
= 0. 



Therefore, we define matrices h^q i"-* = hpqr' ~ h! 



o,(n)' 
xyz—Q ' 

which fulfill the constraint. Substituting h^^V''' into the 

o,(n)" 
xyz—^ 

is obtained. Then the final n-th approximation of the 
direct correlation function is chosen to be 



o,(n)' 

qv il 

from which c" 



a c 



(")" 



+ (1 - a) c>xyz , xyz = , 



,(»)' 

-xyz , 



otherwise . 



(58) 



Note that the self part of the direct correlation func- 
tion can definitely be determined in the form o only (see 
Eq. (|48(l ). In Eq. H58|l . < a < 1 is a mixing parameter. 



(c) \ixyz'' is obtained from \\^qr' by Fourier transform 

and Eq. ^\ yields hi^^^ for xyz 7^ 0. Then, giyz can 

be calculated following App. El Using Eq. (|47|l . ^xyz 

and c^xyz for xyz ^ 0, one obtains new direct correlation 
function matrices, which are mixed to a fraction of a with 

the matrices c^xyz to give ci^t^^ , while c^'^^^j, = c^'^^^^q 
is chosen. 

This prodedure is repeated until a fix point for the ma- 
trices has been reached. Typically, a = 0.1 is chosen to 
avoid divergence. Convergence is assumed if all elements 



of h! 



',(")' 

xyz—Q 



have submerged a certain treshold, which is 



chosen to be 10 times the maximum absolute value of 



any matrix element K, 



o,(n)' 



.AA' 



Additionally, the average 



of all non-zero matrix elements of c'xy^^^ — c^xyz must be 
below a times a second treshold, which is calculated in 
the same manner as the /i-treshold, but by taking also 
the IV = matrix elements into account. It is interest- 
ing that Eq. H48|l cannot be used for the calculation of 
c°^^"2.0' since this always leads to a divergence. 
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After the iteration has converged, also Eq. H48|l can be 
checked to be true. For some systems, only a mixing pa- 
rameter of a = 0.05 or a = 0.02 leads to convergence. For 
some of these critical systems the mixing is not enough 
for the direct correlation function, but also the new total 
correlation function has to be accepted only to a frac- 
tion of a. The medium grey areas in Figs. [51 are the 
areas where the iterative procedure above, for a mixing 
paramter a = 0.02 applied to both the new direct and the 
new total correlation functions, turns from convergent to 
divergent, or begins to yield unphysical negative diagonal 
structure factors (see the next section for explanation). 

Finally, Eq. (|38|l is used to calculate the static struc- 
ture factors. 



C. Numerical results for the correlation functions 

Results obtained from the numerical solution of the 
OZ/PY equations and the MC simulations are presented 
in Figures 0^31 for five different pairs of (a, 6), including 
prolate and oblate ellipsoids. We have restricted the il- 
lustrations of correlators in direct and reciprocal space 
to the matrix elements (/ = /' = 2, m = to' = 0,1,2), 

= 2,/' = 4,TO = to' = 0,1,2) and = T = 4, to = 
to' = 0,1,2,3,4). 

Log-lin representations of the direct space orienta- 
tional correlators Gxyz,\\> are shown along lattice di- 
rections of high symmetry, i.e. {x,y,z) — n (0,0,1), 
{x,y,z) — n (0, 1, 1) and {x,y,z) = n (1,1,1) for n = 

0. 1. . . . , 8 (part a) of Figs. l3ll5|) . Along these directions, 
all Gxyz,\\' are real by the symmetries ()27|l for the cho- 
sen AA'. Note that the PY/OZ results could in principle 
be displayed up to n = 16 and that a step An — 1 corre- 
sponds to different lengths in direct space, namely 1, \/2 
and -s/S for the different lattice directions. In these illus- 
trations, the values n are chosen as abscissas. For each 
m = m! and each lattice direction, a separate picture 
is provided and a logarithmic plotting has been chosen 
for positive and negative values of Gxyz,\\' separately, 

1. e. the negative values are presented as — log \Gxyz,\\'\- 
This plotting shows that the direct space correlations de- 
cay exponentially in most cases. The respective values of 
xyz, AA' are included, too. Note that the scatter of the 
MC data for some AA' and large enough values of n is 
due to the error margins. 

Similarly to direct space, only correlators Sw (q) along 
highly symmetric reciprocal lattice directions are dis- 
played, i.e. (p,g,r) = 5(0,0,7r), (p, g, r) = C(0,7r,7r) 
and {p,q,r) = ^(tt, 7r,7r) for ^ = 0...1, which are the 
correlators from the Brillouin zone center to its edge 
in the respective direction (part b) of Figs. I3I13|I . In 
these figures, ^ is chosen as abscissa and the curves for 
all three reciprocal lattice directions have been put into 
one illustration for each pair AA', which is included in 
each picture. The curves are distinguished by the sym- 
bols A ((0, 0, 7r)-direction), E ((0, tt, 7r)-direction) and A 
((tt, TT, 7r)-direction). All S\\i{q) which are shown are 
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real by the symmetries (I28|l . Additionally, by Eq. H37|l . 
the diagonal elements S'2m2m(q) and S'4m4m(q) are pos- 
itive. The numerically determined correlators Sw (q) 
have been interpolated by cubic splines with the correct 
boundary condition of vanishing gradients for ^ = and 
^ = 1. Note the different scales of the illustrations for 
different AA'. 

Let us start by discussing the results shown in Figs. 13151 
for oblate hard ellipsoids with axes a = 0.4 and b — 1.2, 
corresponding to a packing fraction of (j) ~ 0.3. This 
system is quite close to the MC phase boundary shown 
in Fig. 121 b), but not close enough to find a tendency to 
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4 and m = 



a divergence of some of the Sw (q) . Notice the almost 
perfect agreement of the OZ/PY and MC orientational 
correlators in direct space, in case "where the MC results 
are large enough. In fact, such an agreement appears 
for all investigated oblate ellipsoid systems, for "which 
MC results are available, except for the system "with 
a — 0.72, h = 1.1. Relatively long ranged oscillations 
appear for all correlators along the fourfold lattice direc- 
tion [0, 0, 1] having even m — m! . The other correlators 
along the same direction decay faster and monotonously 
"without oscillation, and also the correlators along the 
other directions decay "without oscillations. Note that 
for a — 0.4 and h — 1.2 the ellipsoids can only interact 
via their nearest neighbours, which are localized along 
the fourfold lattice directions. Therefore, it is tempt- 
ing to assume that the oscillations are primarily related 
to a direct particle interaction via nearest neighbours 
along a certain lattice direction. For denser oblate and 
some prolate systems, the oscillations extend up to many 
lattice constants, indicating the tendency to build up a 
medium range orientational order. For the structure fac- 
tors Sx\i{q), the agreement of OZ/PY and MC results is 
satisfactory. The most significant deviations appear for 
Z = /' = 2 near q = 0. The oscillations exhibited by some 
of the Gxyz,\\' manifest themselves in some maxima at 
the Brillouin zone edge, mainly for correlators S'2m2m(q) 
and 54^,4™ (q) "with m even, but also for the correlators 
S'204o(q)- Increasing a and/or b for oblate systems, the 
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FIG. 6: Same as Figure|3] but for prolate hard ellipsoids with 
axes a = 1.6, b = 0.6. 
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FIG. 7: Same as Figure El but for Z = 2, = 4 and prolate 
hard ellipsoids with axes a — 1.6, b — 0.6. 



OZ/PY results for some of these zone-edge maxima sho"w 
a tendency to diverge, accompanied by a simultaneous di- 
vergence mainly of the remaining I = I' = 2 correlators 
at the zone center. 

Next, "we have chosen the prolate ellipsoid system "with 
axes a = 1.6, b = 0.6, having a packing fraction of 
4> ~ 0.3, for "which some results are sho"wn in Figs. El 
and 13 For this system, oscillations appear additionally 
along the two other lattice directions (i.e. [0, 1, 1] and 
[1,1,1]), which confirm the assumption that they may 
be caused by direct interaction via nearest neighbours for 
appropriate values of b and not too large a. Again, the 
MC results in direct space match the analytical results 
very well, though latter overestimate some of the corre- 
lation lengths. This is clearly seen, for example, from 
the decrease of the correlators Goon,2i2i for increasing n, 
which is different for MC and OZ/PY results. The much 
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FIG. 8: Same as Figure|3| but for prolate hard ellipsoids with 
axes a = 2.0, b = 0.48. 
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FIG. 9: Same as Figure |21 but for / = 2, = 4 and prolate 
hard ellipsoids with axes a = 2.0, h = 0.48. 



too large OZ/PY results for the correlators 6*2222 (q) at 
the zone boundary along the A-direction are perhaps also 
due to this overestimation. OZ/PY pretends an intersec- 
tion of the curves for S'204o(q) along lattice directions A 
and E at ^ « 0.8, which is not present in the MC results. 
Note that the correlators Gonn,2i4i and G'„„„^2242 deviate 
significantly from the regular decrease inherent to other 
correlators. 

Increasing a to 2.0 and decreasing h to 0.48 leads to 
w 0.24 and the results of Figs. I8I10I Now, only one 
of the I = I' ^ 2 direct space correlators still shows 
pure oscillating behavior (Gonn,2222)j and also for the 
other values of 1,1' the number of oscillating direct space 
correlators has decreased. Interestingly, some correla- 
tors oscillate for small values of n and then start to 
decay monotonously (for example Goon,2020, Goon,2222, 
Goon.2242)7 whilc othcrs oscillate from n = 2 on (for 



10" 
10" 
10" 
-10" 
-10" 
-10"' 
10" 
10" 
10" 
-10" 
-10" 
-10" 
10"' 
10" 
10" 
-10" 
-10" 
-10" 
10" 
10" 
10" 
-10" 
-10" 
-10" 
10" 
10^ 
10" 
-10" 
-10" 
-10" 



,j 00n,4444 , 

"', a n 

- 1 

1 . . . 


jj 0nn,4444 [ 

a □ 


1 nnn,4444 

^Og naa 


El 


□ °° ° 


a DO 


(, 00n,4545 1 
; V oa 


] Unn,4i4a [ 




1 nnn,4ji4ji 

1 fl □ □□ 
1 /'^ ,* 

is D ° 


"', S 1? an 

-\.iU\ 

J is a- 


■ ' ' ■ ■ ■ 

^a* Dp 


1 ririri,4242 


ti 00n4141 1 

t J IJ ± -r ± 1 

'\f 


j] 0nn,4141 [ 


1 n nn ,4 14 1 


-I 




D 


^ DO □ 


,j 0Un,4O40 1 
"; I na a 

-vi\ 


J 0nn,4U40 , 
\ V, .DO □ 

A. 


1 nnn,404y 








2 4 6 8 C 
« 


) 2 4 6 8 C 
« 


2 4 6 8 
« 



.40 



4^5 A 



A./'/ 



FIG. 10: Same as Figuregl but for I = I' = 
and prolate hard ellipsoids with axes a — 2 



4, m = 0,1,2,3,4 
,0, 6 = 0.48. 



example Go„„,2020, Go„„,204o)- Though not few of the 
direct space MC correlators differ by sign from the 
OZ/PY correlators (for example Goii,2i2i), the results 
still match satisfactory. Indeed, for some of the corre- 
lators S'4m4m(q) OZ/PY predicts crossings of the curves 
for different lattice directions, which are also present in 
the MC results. Coming back to the direct space orienta- 
tional correlators, more and more of them have become 
non-regular functions of n, and also significant deviations 
from the almost exponential decay, which is present in 
many of the direct space correlators of the other systems 
discussed above, are encountered (see, e.g., Gonri,2i2i)- 
But the MC results for the correlators Goon,4i4i, for ex- 
ample, prove that these non-regularities are not just an 
artefact of the analytical approach. Surprisingly, the 
OZ/PY correlator Go33,4i4i « -2.97 x 10"^ is to small 
to be displayed in the pictures, though the surrounding 
correlators (i.e. GonTi,4i4i, n — 0,li2,4,5,6) clearly are 
not. 

Now we pass to more and more elongated prolate el- 
lipsoids and investigate what happens. In the limit of 
this process one obtains the hard-needle system, which 
is discussed, e.g., in 18] . Ellipsoids with axes a = 3.6 
and b = 0.24 have a huge aspect ratio of = f = 15. 
Of course, the packing fraction of these ellipsoids on the 
sc lattice is quite low, about 10%, but the frustration 
effect of the rigid lattice may provide completely new ef- 
fects in comparison to a liquid. As shown in Fig. 1111 the 
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FIG. 11: Same as Figure |3 but for prolate hard ellipsoids 
with axes a — 3.6, b — 0.24. 
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behavior of the correlators Gxyz,2m2m is more regular in 
comparison to a = 2.0 and b = 0.48. Oscillations have 
disappeared completely. The OZ/PY results clearly un- 
derestimate the S'2m2m(q)-correlators for small q. This 
is also the case for the 5*2^2771 (q)^correlators of all other 
investigated systems in the neighbourhood of a = 3.6, 
h = 0.24. 

The last system "we present consists of ellipsoids of 
axes a ~ 4.8 and h = 0.24, for -which it is Xq — 20 
and (j) w 16%. The results of Figs, d and US belong to 
this pair (a, 6). The correlators in direct space essentially 
sho"w up monotonous decay, and the correlation lengths 
have clearly increased in comparison to the elliposids of 
axes a = 3.6 and h = 0.24. Unfortunately, this system 
lies beyond the MC phase boundary (see Fig.l^L), so that 
no MC results are available. Despite of the monotonous 
behavior of most correlators, the Goon,2020"Correlators, 
for example, still sho"w non-regular behavior, "which van- 



ishes for the most part S"witching to a = 5.6 and h = 0.24 
(these results are not shown here). In Fig. ^1 the be- 
ginning of a divergence of the S'2i2i(q — 0)-correlator is 
seen, the corresponding value of the (a — 5.6, h — 0.24)- 
system being « 20. Also, "what can be seen from Fig. 1131 
the / = 2, r = 4 correlators have got significant struc- 
ture, and they seem not to diverge like the I = V = 2 
correlators. 



V. DISCUSSION AND CONCLUSIONS 

Our main goal has been the study of static orienta- 
tional correlation functions for a molecular crystal in its 
disordered phase. For this, "we have derived the OZ equa- 
tion, ■well-kno"wn in liquid theory, for a rigid periodic lat- 
tice with internal orientational degrees of freedom. As a 
closure relation we have adopted the PY approximation. 
As pointed out, there are differences of the present ap- 
proach to that for liquids. One of them is the fact that the 
OZ equation only involves the direct and total correlation 
functions c°„,(f7,f}') and /i°„,(f^, f7'), whereas the PY 
approximation relates c„„'(f2,r2') and /i„„/(ri, SI'). The 
functions with the superscript o do not contain a constant 
part with respect to f7 or fi'. Another important differ- 
ence is the one-particle orientational distribution func- 
tion p(^^(ri). In the isotropic phase of a molecular liquid 
it is p(^^(ri) = but due to the anisotropy of a crystal 
p^^-'(ri) exhibits a nontrivial r2-dependence. In order to 
solve the OZ/PY equations, one has to calculate p'^'^\^) 
separately. In our case, we have performed MC simula- 
tions. Analytical approaches are also possible, e.g. for 
fixed a one could perform a kind of virial expansion for 
small h. 

Despite these differences, the form of the OZ equa- 
tion for a molecular crystal is quite similar to that for 
molecular liquids In order to explore the applicabil- 
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ity of the lattice OZ equation in combination with the 
PY approximation, we have solved these equations for 
hard ellipsoids of revolution on a simple cubic lattice. 
Due to the orientational degrees of freedom, the orien- 
tational correlators Gnn',\\' in direct space or S'aa'(i) 
in reciprocal space, with A = (/to), are tensorial quanti- 
ties. Accordingly, the self consistent numerical solution 
of the OZ/PY equations requires a truncation at /max- 
We mainly have chosen /max = 4. As a result, we have 
found orientational correlators which have less structure 
in direct and reciprocal space than for liquid systems. 
Nevertheless, there are some interesting features depend- 
ing on the length a and b of the ellipsoid axes. For oblate 
ellipsoids and prolate ones of large enough 6, some of the 
direct space orientational correlators exhibit oscillations 
in certain lattice directions. Since the oscillations have 
period two they lead to maxima of S'aa'Cq) at the Bril- 
louin zone edge for some AA'. Although no long range 
orientational order exists, the oscillatory behavior origi- 
nates from an alternating reordering of the ellipsoids on 
a finite length scale, which can extend up to many lattice 
constants. 

Decreasing for prolate ellipsoids b and increasing a 
leads to a disappearance of almost all of these signifi- 
cant oscillations. In this case, the correlators S'2m2m(q) 
take their absolute maxima at the Brillouin zone cen- 
ter (cf. Figs. [TT| and I12|l. indicating nematic-likc fluc- 
tuations, while the same behavior is not found to the 
same extent for the other correlators (i.e. S'2m4m(q) (cf. 
Fig. El and 54^4™ (q) (cf. Fig. nHJ). The behavior of 
the correlators 5'2m2m(q) resembles a liquid of ellipsoids 
with aspect ratio larger than about two which forms a 
nematic phase fl9j . Surprisingly, increasing a for fixed 
h more and more the OZ/PY results for S'2m2m(q) lead 



to a divergence at q = 0, which indicates the tendency 
to establish a long range nematic-like order. This find- 
ing demonstrates that PY theory can yield the onset of 
a phase transition to an ordered phase as it was already 
found before for a liquid of hard ellipsoids 4'j . 

Some correlators for appropriately long prolate ellip- 
soids show up highly irregular behavior in direct space, 
which is nevertheless consistent with MC results (where 
available) and therefore has to be taken serious. We sup- 
pose the frustration effect of the lattice to be the reason 
for this non-regular behavior. 

Comparison of the PY results with those from MC sim- 
ulations shows a satisfactory agreement. But the quality 
of this agreement is less good as it is, e.g. for a liquid of 
hard spheres. The reason for this may lie in the PY ap- 
proximation. For a liquid its physical content has been 
elucidated by Percus (see also ref. 0) by use of a 
grand canonical ensemble. Since in a molecular crystal 
the particles are fixed, it is not obvious how this rea- 
soning can be used. Of course, it might be interesting to 
investigate whether other closure relations [Jl or the turn 
to much higher /,„ax can lead to an improvement. 

To conclude, we have demonstrated that an extension 
of the OZ equation in combination with the PY approxi- 
mation to molecular crystals leads to satisfactory results 
compared with MC data. 
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APPENDIX A: CALCULATION OF /i^„,(J7,n') FROM h'^;°{Q.,Q.'), n ^ n' 

In Apps. I AID I we discuss the most general case K C S^. U K = S^, the superscripts K which are used have to be 
dropped to be consistent with the rest of the paper. 
Using definiton (|5()(l . we have with Eq. (|51a|l 

h^„,inM)^h^;^%n,n') + j^Jh^„,in,n')dn + j^Jht,in,n')dn'-j^ JJ h^„,in,n')dndn' (ai) 

K K KK 

If we take //^ p'^^ i^) ■ ■ ■ p'^' {^') dfl' on both sides of Eqs. (jXTJ, it remains due to Eqs. (O 



IK 

K,o 



KK KK 



h^„, {n,n') dn dn' = / / p^^) (n) /i^„? {n, n') p^^'> (n') dn dn' . (A2) 



Only . . . p^^^ i^') dQ' apphed to both sides of Eq. (|XT)l yields 

^ ^ i^^^')dn' ff i^, ^')dn dn' - I /if„? {n, n') p^) (n') dn' (A3) 



K KK K 



(analogously, p'^^^n) . . . dn must be applied in the same manner). Then, ljA2p can be used with the first term on 
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the rhs of ljA3p . and finally we find that 



K 



K 



KK 



KK 



dn" dn" 



(A4) 



Therefore, /i^^, (fi, ft') is determined by h^^{rt, ft') via Eq. (|A4|) . It is easy to prove that this h^^,{n, ft') is the only 
one leading to the given and simultaneously fulfilling Eqs. I|15|l . i.e. /i°'^(ri, f2') determines h^^,{n,fl') 
uniquely. 



APPENDIX B: CALCULATION OF THE 
MATRIX ELEMENTS g^^.^xx', ^ ^ n' 

For the PY approximation in matrix form 147() . the 
matrix elements .9,^,,^;^, = lx\'+f^nn',\y are needed. In 
the following, it is convenient to introduce the projection 
operator = [P'^f = {P^)\ defined on 5^ x 5^ as 

p-i^M) = \m^').ji.^:^K., 

^ ' 1 , otherwise 



P 



K 
W 



(Bla) 

Y^{VL)Yy{VL)dQ., (Bib) 



K 



and 

d{n,^') ^ At: p'^^\n)5{n\vt') 



(B2a) 



dxx' = 47r i 



I'-i / Ji) 



p^'>{n)Yx{n)Yx'{n)dVL. (B2b) 



K 



Using (IbU, the A transform of D{Q.,^') (see Eq. (tT^ ) 
can be represented as 



Next, Eq. (|A4|I is used. The matrices h^ ? are deliv- 
ered by the OZ equation. According to l|B3(l . the matrix 
belonging to //^ (f]) h^^° {n, n') (f7') dO. dn' is 



(r!) ;if„? (r!, r!') (r!') dn dn' 



CK 



— A-TT 

— ^~ ^A.OO ^00, A' 



KK 



(n) /i^„?(r!, n') (r!') rfr* dn' , 

(B4) 



where by use of ljB2b|l one easily computes that (using 
the summation rule) 



KK 



(n) /if„? in,n') p'-^^ {n') dn dn' 



(B5) 



D 



XX' — "AA' 



d 



A, 00 "00, A' 



(B2c) 



The matrix elements (/^ h^;^°{n, ft') p^^Hfl') df7')AA' are 
also straightforward. Since the angular function depends 
not on Q' , ^/in P^^ y must occur, as in Eq. (|B3(I . The 
Though defined on x S"^, d = d^ can be assumed, dependence on f2 and the integrals are represented by 
since d{n, n') = on {S^ x S^) \ {K x K). The inverse 

fimction for d on is ''pW(n) ' ^^te that for K ^ 

one has P^y = Sx,oo Soo,X' ■ 

Now we calculate ^y. First, the matrix elements 
1^^, are given 



4"" ^nri/ XX" "^a'^oOj SO that we finally have 



hl';;^°{n,n')p^'\n')dn' 



XX' 



1^ 

-^AA' 



1 1 r; {n) Yx' {n') dn dn' = An p^^,^ p^l^' ■ 

(B3) 



- h^'° rl P^ 

~ "nn',AA" "A" ,00 -Hio, A' ■ 

(B6) 



KK 



Altogether, the matrix elements h^^, ^y, which must be 



completed by (jB3p to give g. 



.K 



ynn' ,XX' ' 



nn' , AA' 

are 



hnn',XX' - ^nn'.XX' ~ ^nn' .XX" ^A",00 ^TO,A' ^ -^^00 C^OO.A" h^n' ,X" X' + ^A^OO C^OO.A" h^n',X"X"' ^A"',00 Pm,X' 

= {Pxx" - Px,oodoo,x") ht\x"X'" {Py-x' - Polx'dx'" ,00) (B7) 



18 



APPENDIX C: CALCULATION OF f^„',\x', w / n' 



Contrary to the lattice correlation functions, fnn'{^, is not affected by the lattice and refers exclusively to what 
happens between two particles. Therefore, it is advantageous to use the r-frame for the calculation of the matrix 
elements, i.e. the coordinate system in which the connecting line of the particles coincides with the z-axis, and later 
to transform them back to the laboratory system. For K — in the r-frame, the matrix elements have similar 
properties as in p^: 

j^r A ( 1 V^"^ f^* {C^^\ 

Jnn' ,lm,l'm' J nn' ,11' m '^mm' i Jnn'M'm I ^) Jnn',ll'm : Jnn',ll'm J nn' ,IV —m ■ V^-*-/ 

In the following, the abbreviation Qim{6) — (— 1)™(^^)^ ( (l+mjl )^ sin 6> (cos 6>) is used, where we define 
Pi.-m{cos9) = (— 1)™ PimicosO) to covcr all possible values of m. Then it is 

fL'll'm=^''~' [ I I f r{\^nn'\,0r,9'^,cl,'^-cl,r)Qlrn{0r)Ql'rn{0'r)e''^^^'-^-^^ (C2) 



lnn',ll'm — '- J J J J W-^nn' \, " r, " r-^ Vr ~ Vr I ^ lm\^r ) ^ I' m\'J r ) ^ niy^ "i^r "V^r "V'r 



Switching to the new variables = — (j)r and (/)2 = 0^ + 'Pr brings about a functional determinant factor of ^ : 

4ir 'I'(02) TT IT 

-*(02) 

where *((/>2) = 02, < 02 < Stt, *((?!)2) = 47r - 02, 27r < 02 < 47r. Now, the symmetries /''(|x„„/|, 6*^, 0i) 
f''{\xnn'\,Or,0'j.,-(l)i) = f \ , r , 0'^ , 2tt - 0i) Can bc uscd: 

in *(02) IT TT 

fnn'M'm=^'''' J J J J ri\^nn'\,0r,9',,cj)l)Ql„^i9r)Ql',ni9r) COsim(j)i)d9rde', d^ld^2 




® 

2ir <p2 27T 2-K 2-K 2-K 

= 2i''"' y" J ®d<j)i d(j)2 = 2i''"'y" j ® d02d0i = 2z''"' ^ (27r - 0i) ® d0i = 27ri''"' ^ ® d0i 

01 

TT TT TT 

= A7:i''''J J J r{\^r.n'lOr,0'r,MQlrn{Or)Ql'rn{0r) COs{m4,i)d9rd9',d(lji. (C4) 


For particles of inversion symmetry, only the (II' even) matrix elements do not vanish, and by using 
f'^{\xnn'\,0r,9'^,<j)i) = /''(|x„„/ 1 , 6*^, TT - 6*^ , TT + 0i ) = . . . for that case we find that 

'2 

/™',«'™ = 167ri''-'y J J ri\^nn'l0r,9',,^i)QiU0r)Ql'mi0'r) cos{m^,)d9rd9'^d(l),. (C5) 


Finally, having transformed back to the laboratory frame and the matrix elements fnn',\x' in hand, they can be 
projected due to f,^,, = f„„- P^. 



APPENDIX D: CALCULATION OF c/^'^,{n,n'), 
n^n', FOR HARD PARTICLES 

The appendices IdI and IeI are given in terms of angular 
functions only. The discussion in terms of matrices is 
much more complicated. 



For hard particles, the values of 0^^/(17,0') in the ar- 
eas of overlap remain undetermined by Eq. (|45|1 . and the 
question arises if these values are also uniquely deter- 
mined by c^jj? (il, f2'). That this indeed is true will be 
proven below. 

Most important for the following considerations is that 
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no space angle fif, S if can be found for which the pair 
(nfi, uTIq) has overlap for all U G K. Otherwise, it would 
be J7g ^ K. Note that, contrary to the case of overlap, 
for appropriate size of the particles and special values of 
JIq G K, the pair {nil, tiTIq) has no overlap for all fl E 
K. Given a solution c^^, for which c^^?(f7,17') = 
(ii'^cf„,ii'^)(r!,l)'), it is cl,(r!,r!') = 0, if [n^^n'^') 
has no overlap. Eqs. (|50|l . H51a(l and the symmetries of 
the correlation functions given in Sec. Ill Bi show that 



V 



K.i 



n'n 



in') 



(Dl) 



Let us fix Vt' ^Vl'^eK and let K(Vt'^) d K he the area 
such that {nVt^n'Q!]) has no overlap for all Vl E K{il'i). 
Then it follows from Eq. l|Dip 



(D2) 



for all n E K{n[), because of c^„,(f^, , f^i) = 0. Since 
c^^?(ri, ri'^) is known, we have found c^„, (fi) on K{n[) 
up to a constant depending on n[ . Next we choose ft' = 
E K and fl'^ ^ n[. Then c,f„,(fJ) is obtained for 
E K{fl'2) up to a constant depending on fij. This can 
be continued by choosing fi^, i = 1,2, ... ,M such that 
LiiK{V,[) — K. In that case one can construct c^, (fi) 
on K, for example in the form c°:^, (£7) having vanishing 
constant part. The same is done for c°',^(r2'). Finally, 
we have 



(D3) 



where c^„, follows immediately if ljD3|) is used with any 
Q,,Vt' such that {nVl,n'U,') has no overlap. 



lim L>^(f],17') = L>°°(17,rj') 



(E2a) 
(E2b) 



This is true for since the one-particle den- 

sity varies continuously with V. If also all 

other functions exist in the limit = 00 and if the 
limit interchanges with integration, we can write down 
Eqs. H43a|l . iflljl . (|A4|I and H32a|l including the func- 
tions hlA^,n'), hl:,{n,n'), cX„,in,n'), C(r!,i7'), 

D^{n, n') and (f], rj') for F < cx) and take this limit. 
Then it is easy to see that the limiting direct and total 
correlation functions are a set of functions /i^, (il, £7'), 
c5^,(ri, ri'), i.e. they solve the OZ/PY equations for 
F = 00 on 52 X 52. Note that, as hX„,{n, O'), h^„,{n, n') 
fulfills the constraints given by Eqs. p5|) . 

Now we investigate if a solution /ij^, (fi, f2') and 
c^, (ri, ri') leads to the correct static structure fac- 
tors Eqs. P7|l . for F = 00. For this, since 

{R^ £)oo f]')^ with R^{n, n') due to Eq. we 

have to prove that (i?^ i?^)(f7, f}') = /i^„?(f7, f}'), 
n, fl' E K , where the former symbol o refers to removed 
constant parts with respect to and the latter with 
respect to K . This can be done as follows. 

a) It is tempting to assume that the functions 
h'^„,{n,n') and 0(^2,^1') restricted to K mat ch the 
functions /if„, [n, n') and cf„, {n, n') of Sec. mTRl Then, 
the PY approximation H44|l restricted to K is fulfilled 
trivially. Additionaly, the restricted function ftj^, (il, Q!) 
fulfills the constraints of Eqs. p5|l (rewritten for K 
instead of S*^), since these integrals do not change if 
one restricts the integration domains to K (because of 
p'^^\n) = on k). 

b) Because of 



APPENDIX E: SOME REMARKS ON HOW THE 
CASE K CS'^ CAN BE TREATED 

In this appendix, we want to answer two questions con- 
cerning the case of hard core particles, whose geometry 
leads to K G . First, assuming finite interaction if the 
"hard cores" of such particles overlap, i.e. 



Vnn' = < 00 (overlap) , 



(El) 



what happens to the solutions {Vt, D!) and c^„, (57, U!) 
of th e OZ /PY equations l|4aa|l and igU, together with 
Eqs. (|A4|) and (|32a|l . which are defined on 5^ x S*^, if we 
take the limit V — > 00 ? Second, does a S*^ x S'^-solution 
/i^„, (17,17') and C„,(17,17') of the OZ/PY equations for 
V — 00 lead to the same static structure factors as the 
solution /if„, (17, 17') and c^„, (17, 17') of Sec. ITITBt ' 

The first question can be answered using the fact that 
the function 17^(17,12') (cf. Eq. iCHIl) and the Mayer 
function Eq. 145|) behave continuously if we take the limit 



j (^5(17|17') - I y" (^<5(17'|17") - ^ /(17") dl7"| dl7' 

(E3) 



(5(i7|i7') - — ) f{n')dn' 



K 



is the removal of the constant parts of a function with 
respect to S"^ and a subsequent restriction to K and re- 
moval of the constant parts with respect to K equivalent 
to restricting the function to K first and then to remove 
the constant parts with respect to K. 

c) Next we notice that the OZ equation H43a() for x 
decouples for V — 00 ox\ K y. K completely from K, 

as does the OZ equation of Sec. IIII 51 

d) If we write down this decoupled OZ equation with its 
solution /i^',°(17, 17') and c^°(17, 17'), restrict it to x if 
and apply the projection operator i?^(17, 17') from both 
sides, it involves the functions [R^ K^n° ^') ^^^d 
(i?^ 0° R^ W, ^')- Due to b), this is the OZ equation 
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on K X K, written with the functions fi') and 

c^,(r2, ri'), which are restricted to K and have removed 
constant parts with respect to K. 

Now we have shown that the initial assumption made 
in a) was right. Therefore, to solve the OZ/PY equations 
for the case K C S^, one can solve these equations also 



on 5*^ for V = oo, since the irrelevant K'-parts of the 
involved functions are projected out. We have also made 
plausible that the limiting process from finite to hard 
core interaction leads to such a solution on S*^, and that 
therefore the structure factors behave continuously for 
V — *■ oo. This is exactly what one expects by intuition. 
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